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ABSTRACT 

We present a multi-wavelength analysis of the history of star formation in the 
W3 complex. Using deep, near-infrared ground-based images, combined with im¬ 
ages obtained with Spitzer and Chandra observatories, we identihed and classi- 
hed young embedded sources. We identihed the principal clusters in the complex, 
and determined their structure and extension. We constructed extinction-limited 
samples for hve principal clusters, and constructed K-band luminosity functions 
(KLF) that we compare with those of artihcial clusters with varying ages. This 
analysis provided mean ages and possible age spreads for the clusters. We found 
that IC 1795, the centermost cluster of the complex, still hosts a large fraction 
of young sources with circumstellar disks. This indicates that star formation was 
active in IC 1795 as recently as 2 Myr ago, simultaneous to the star forming ac¬ 
tivity in the hanking embedded clusters, W3-Main and W3(OH). A comparison 
with carbon monoxide emission maps indicates strong velocity gradients in the 
gas clumps hosting W3-Main and W3(OH) and show small receding clumps of gas 
at IC 1795, suggestive of rapid gas removal (faster than the T Tauri timescale) in 
the cluster forming regions. We discuss one possible scenario for the progression 
of cluster formation in the W3 complex. We propose that early processes of gas 
collapse in the main structure of the complex could have dehned the progres¬ 
sion of cluster formation across the complex with relatively small age diherences 
from one group to another. However, triggering ehects could act as catalysts 
for enhanced efficiency of formation at a local level, in agreement with previous 
studies. 
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1. Introduction 


In the current picture of star formation, Giant Molecular Clouds (GMC) are highly 
inefficient factories, in which only a small fraction of the available gas is converted into stars. 
Moreover, star forming regions in GMC complexes are highly heterogeneous. There is a 
considerable diversity among the stellar aggregations they produce. Differences are evident 
even at the level of their basic properties, e.g. sizes, numbers and density structures. 


Aggregations of young stars in molecular clouds are classihed as Embedded Star Clus¬ 
ters. These clusters are born as bound systems of gas and stars. Once gas is removed, 
embedded clusters may or may not stay bound. They either survive (e.g. open clusters), 
or dissolve to become part of the held population. For a majority of clusters in the galaxy. 


the l atter is the common outcome (referred as the “infant mortality effect”, iLada fc Lada 
2003h. However, recent discussions consider that cluster morphology is more heterogeneous 
flBressert et al.ll2010l: Ide GriisI 1201111 . Num erical studies suggest that small groups of stars 


may be abundant flAdams fc MversI 120011) and they may eithe r merge to forn i clust ers or 
disrupt rapidly into the held. According to studies like that of iKruiissen et ahl (120121) . star 
clusters may assemble from sub-clusters that merge into larger entities after the gas from 
which they formed is both disrupted by stellar feedback and torn apart by tidal shocks from 
the surrounding cloud (the so called “cruel cradle” ehect). 

Unfortunately, the dynamical picture of cluster evolution provided by numerical studies 
is difficult to directly compare with observations. On one hand, the initial conditions of 
cluster formation, a key requirement for realistic simulations, are still under debate. On the 
other hand, young cluster evolution from observations, would require of additional informa¬ 
tion on the kinematics of stars (e.g. radial velocities), and spectroscopic age estimations 
of individual sources to infer the professio n of formation. Suc h studies are only beginning 
to arise from surveys like APOGEEpJ (e.g. iGottaar et ahl 120141) . and are limited to nearby 


{d < 1 kpc) regions. Meanwhile, it is possible to partially reconstruct the history of star for¬ 
mation in a region from photometric information, which can provide evolutive classihcation 
and the spatial distribution of young sources, and from molecular gas emission maps, which 
can provide gas distribution and kinematics. 

Surveys of GMGs show that embedded star clusters are rarely (if ever) born in isolation. 


^ Sloan Digital Sky Survey, Apache Point Galactic Evolution Experiment 
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Most embedded clusters form as part of “families” related to a particular complex, defining 
together a history of star formation from various levels of interaction. For instance, massive 
stars of one cluster may have influence on the efficiency of star formation of a neighboring 
cluster in the same cloud. Considering all this, it should be clear that star formation is 
highly dynamic; the embedded populations we observe are snapshots of a rather convoluted 
process of evolution and interaction, that changes signihcantly along the star forming history 
of a complex. 


In order to investigate the progression of cluster formation under the influence of the 
local environment, we considered W3, a prominent cluster forming region in the “Heart and 
Soul” molecular complex. W3 is considered a clear example of sequential forma tion, where 
clust e r formation was indu ced by the expansion of the giant HII region W4 flLada et ah 


19781: iThronson et al.lll985l) . One of our goals is to trace the formation of distinct embedded 
clusters in this region and to attempt to reconstruct the star forming history of the cloud. 
The W3 Complex is located at a distance of 2.04T0.107 kpc B, and has hosted at least 
three major episodes of recent cl uster formation. According to spectroscopic studies, IC 


1795 formed hrst, about 3-5 Myr fjOev et al.ll2005l) ago, followed by the W3 “main” cluster 
locat ed to its West edge, and the W3(OH ) cluster group to the East, both with ages of 2-3 
Myr flBik et al.ll2012l:lNavarete et al.ll20111) . It has been suggested that IC 1795 triggered the 
other two episodes in a hier archical progression flOev et al.ll2005l) . The Chandra study by 
Feigelson fc TownslevI (120081 ) suggested that the clusters in W3 extend widely and are highly 
structured, with sources located at relatively large distances from the dynamical centers, 
including a relatively isolated O star that might have escaped from the main cluster. 


In two recent studies, iRivera-Ingraham et al.l (120111. l2013l) made use of mid-infrared 
photometry and far-IR emission mapping from the Spitzer and Herschej^ space observatories, 
that comprise the entire W3 region. They were able to catalog hundreds of young stellar 
sources (YSOs) across the complex and determine their spatial distribution. They compiled 
important evidence that projected distances among YSOs are consistent with cluster forming 
clump scales, favoring the cluster forming mode. They also concluded that small aggregations 
and distributed populations, account for a signihcant fraction of the recent stellar production 
in the region. They suggest relative large age spreads from central to external regions in 
W3 and proposed a “convergent constructive feedback” scenario, where the gas flows from 
massive star formation clumps and collects into new dense regions, favoring a progression of 


^based on Very L ong Baseline Interferometry (VLBI) parallaxes to H 2 O masers in W3(OH) by 
Hachisnka et al. ( 200(lll 


^Herschel is an ESA space observatory with science instruments provided by European-led Principal 
Investigator consortia and with important participation from NASA. 
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Fig. 1.— A false color mosaic made from J (blue), H (green) and K (red) images from our 
OMEGA 2000 observations. . 
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formation and the age spread. 

In this study we emphasize the relative importance of different young stellar groupings 
in and around IC 1795, W3-Main and W3(OH), augmenting the level of detail achieved in 
previous studies. We present a new set of deep, high resolution near-IR photometry that 
we combine with other available datasets, allowing us to increase the number of young star 
candidates in the complex. We use a K-band luminosity function analysis to investigate the 
presence of an age spread, and to help reconstruct the history of star formation in the cloud. 
Finally, we investigate gas kinematics near the clusters, providing additional information on 
the gas-star interaction during the early evolution of the complex. 


2. Observations and Data Reduction 
2.1. Near-infrared imaging 

We obtained near-infrared images of four helds in W3, covering a large area (0.58 deg xO.46 deg). 
Images were obtained with the Omega 2000 camera at the 3.5m telescope in the Calar Alto 
Observatory of the Centro Astronomico Hispano Aleman (CAHA) in Almeria, Spain. Omega 
2000 provides a 15' x 15' held of view (FOV). Observations were made in J, H and K (1.1, 

1.6 and 2.2 /mi, respectively). A list of all helds observed is listed Tabled which gives the 
held identihcation, the center of held positions, observation date, hlter, seeing (estimated 
from sigma-clipped average full width half-maximum (FWHM) of the stars in each held), 
and the peak values for the brightness distribution, which is a good estimate of the sensitiv¬ 
ity limits achieved. The brightness distribution peaks are in all cases at or above J = 20.5 
mag, H = 19.75 mag and K = 19.25 mag. 


2.1.1. Image reduction 


The Omega 2000 images were reduced with modihed versions of the FLAMINGOS near- 
infrared reduction and photometry/astrometry pipelines, which are built in th e standard 
IRAF Command Language environment. One pipeline fsee lRoman-Zi]nigall2006l) processes 
all raw frames by subtracting darks and dividing by hat helds, improving signal to noise ra¬ 
tios by means of a two pass sky subtraction method, and combining reduced frames with an 
optimized centroid ohset calculation. We used dark frames and dome hats obtained within 
48 hours of each obser vation. The h nal combined product images were then analyzed with 
a second pipeline, (see Levine 2006h. which identihes all p ossible sources from a given held 
using the SExtractor algorithm flBertin fc ArnoutsI Il9961) . We improved the SExtractor 
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detection efficiency in the near-IR images by u sing a Gaussian convolution filter and max¬ 
imum deble nding fsee [Bertin fc Arnouts 1996). The pipeline then performs Daophot PSF 
photometry fjStetsonlll987l) . calibrates observed magnitudes to a zero point and hnds accurate 
astrometric solutions. 

After astrometric solutions were found, individual held images were combined into mo¬ 
saics using Montage 0. In Figured] we show a RGB panorama constructed from a combina¬ 
tion of the J, H and K mosaics. The image represents a complete spatial coverage of the IG 
1795, W3-Main and W3-OH clusters. In the near-infrared the region is transparent to most 
of the prominent reg ions of nebulosity and obscuration observable in optical images (e.g. 
Qgura fc Ishidal 119761) . The near-IR images at W3(OH) are particularly interesting because 
they reveal, with unprecedented detail, several small, embedded stellar groups that lie in 
a small “chain” structur e north and east of the W3(OH) cluster (these g roups have been 
previously identihed; e.g. iFeigelson fc TownslevI 120081: iNavarete et al.ll201ll) . In Figure [2] we 
show a close-up image of this region, where we now combine H and K images with the Spitzer 
IRAG 3.6 yum image to enhance illuminated nebulosity features and highlight the most red¬ 
dened sources. We have labeled the W3(OH) cluster with the letter ‘A’. The two other 
most conspicuous groups in the ’chain’, both associated with B-type stars (INavarete et ah 
201ll) have been labeled with letters ‘B’ and ‘G’. The former group is not as prominent as 
an over-density but as it can be seen in the mosaic, it is associated with thick nebulosity, so 
it is possibly more deeply embedded. There is one more bright source sitting at the center 
of a cavity in between groups B and G, which we suspect could be an additional sub-group. 


2.1.2. Catalog preparation 


The OMEGA 2000 photometry was calibrated relative to the Two-Micron All Sky Sur¬ 
vey (2MASS), using catalogs retrieved from the All Sky Release Point Source databases. 
Final photometry catalogs were prepared with aid of TOPCAT-STIL flTavloiil2005l) . We com¬ 
bined individual photometry catalogs from each frame and hlter and merged catalogs from 
all four helds into a master photometry list. In the overlapping areas, we selected the dupli¬ 
cates with the smaller total photometric error across the three bands in either held. Finally, 
we added 2MASS entries to replace failed measurements from saturated sources. Our hnal 
near-IR catalog contains a total of 72168 sources. 


"‘http: //montage.ipac.caltech.edu 
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Fig. 2.— RGB color mosaic of the W3(OH) region, from Spitzer IRAC 3.6 /im (red) and 
K and H OMEGA 2000 images (green and bine, respectively) showing possible snb-clnster 
strnctnre in W3(OH). The W3(OH) cluster is labeled as “A”, while two possible sub-clusters 
are labeled “B” and “G”. These last two labels are placed closer to B-type sources. 



2.2. Mid-infrared Spitzer imaging 


We performed a new reduction and photometric extraction of the W3 mosaics from 
project 1-127 available in the Spitzer Heritage Archive. Mosaics were constructed by pro¬ 
cessing both the short and long exposure (0.4 and 10.4 sec) Basic Calibrated Sets (BCD) 
from IRAC (all four bands) and the BCD set from MIPS (24 /rm only) observations using 
the MOPEX packag^. We used the IRAC short exposure datasets to improve extraction 
of sources down the centers of W3-Main and W3(OH) regions, which present large areas of 
saturation (increasing with wavelength) due to bright nebulosities in the long exposure sets. 
However, even at the short exposure images, the bright nebulosity made it very difficult 
to obtain uniform quality from aperture photometry in the IRAC images, so we performed 
PSF photometry with Daophot on detection lists obtained with SExtractor. In this case, we 
improved the d etection efficiency by us ing a “Mexican hat” convolution hlter and reduced 


deblending fsee iBertin fc Arnoutslll996l) (compare to the near-IR, in section [2.1.11 above). 


Then, our pipeline selects a group of moderately bright and relatively isolated sources in 
order to make a PSF model for each Daophot run. In all cases we were able to £t the PSF 
model to a large majority of the detected sources. Final source detection lists were slightly 
cleaned with the aid of Phot Vis flGutermuth et al.ll2004l) to remove spurious detections. 


The resultant photometry lists showed a net reduction of up to 12% in the photometric 
scatter to the aperture photometry obtained with the APEX software. For the MIPS 24 
/rm mosaic we used the APEX aperture photometry pipeline, which was good enough for 
non-saturated regions of the field. In heavily saturated regions, particularly near the centers 
of W3-Main and W3(OH), missing pixels impeded the ability to make good measurements, 
and we could not improve these measurements with PSF photometry. Zero points for our 
Spitzer photometry catalogs were checked by direct cornparison of our photometry lists with 


YSO candidates in 

Rivera- 

neraham et ah 

(2011 

). 

and 

Rivera-Ineraham et ah 

(2011 

) is excellent and 


The hnal photometry catalogs from the IRAC and MIPS 24 pm mosaics were then combined 
into a mid-IR catalog with 37267 sources, which we later merged with the near-IR catalog 
from our CAHA dataset to form a master source catalog. 


® http://irsa.ipac.caltech.edu/data/SPITZER/docs/ 
dataanalysistools/tools/mopex / 

®only IRAC [3.6] and [4.5] bands are available as this is a warm mission extension of the original Galactic 
Legacy Infrared Midplane Survey Extraordinaire (GLIMPSE) project 
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2.3. Other datasets 

In this study we also make use of other datasets available in literature and public, 
web-based archives: 


1. From the Chandra Source Catalog (CSC; Evans et ahl 120101) the positions and X-ray 
photometry of a total of 611 sources in a box of 30' around the center of IC 179^. 


2. The ^^ CO(2-l') and ^^CO(2-l) maps of the W3 region from the study of iBieging fc Peters 

feniih . 


3. The 2.5-level (science grade) SPIRE and PACS mosaics of the W3 region from the 
Herschel Space Telescope data archive. 

4. The B qlocam 1.1 mm mos aic of W3 from the the Canadian Galactic Plane Survey 
(GPS; Aguirre et ah 2011 ). 


3. Data Analysis and Results 
3.1. Identification of Young Stellar Sources 


Using our master infrared photometry catalog we identified Young Stellar Objects 
(YSOs) in the W3 complex. First , we identified Class I/O and Class II sources from IRAC 
colors using the criteria applied by Ybarra et ahl (120131) . which are in turn based in the IRAC 
color criteria of Gutermuth et ah ( 20081) . Then, using the Chandra CSC catalogs we identi¬ 
fied Class III sources as X-ray sources with infrared counterparts and no excess associated 
to circumstellar material. Finally, we identified additional Class I and Class II sources as 
X-ray CSC sources with infrar ed excess, by applying c riteria that combine near-infrared and 
IRAC colors (see appendix of iGutermuth et ahl 120081) . 


Figure |3] shows a. H — K vs. K — [4.5] color-color diagram for all infrared sources that 
coincide in position with X-ray GSG sources. The colors of the symbols are indicative of the 
Hardness ratio, calculated from the hard and soft X-ray fluxes, H and S', estimated from 
aperture photometry measurements, as {H — S)/{H + S). It is clear that most of the sources 
have large hardness ratios, which indicates that in most cases X-ray sources are obscured by 
large amounts of dust which reduces the soft X-ray emission. 


^Observation ID numbers for W3 ACIS-I data are 446, 611, 7356, 5889, 5890, 5891, 6335 and 6348 
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Fig. 3.— H — K vs. K — [4.5] color-color diagram for CSC X-ray sources coinciding with 
an source in our infrared master catalog. Open diamond symbols are sources without an 
infrared excess indicative of a disk. Open square symbols indicate sources classihed as Class 
1. Open triangle symbols indicate sources classihed as Class II. The colors of the symbols are 
coded to the hardness-to-soft X-ray ratio, as indicated in the color bar. 
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We identified a total of 1008 YSOs in W3, distributed as follows: 60 Class I/O candidate 
sources, 780 Class II candidate sources and 167 Class III candidate so urces. Our list contains 


a sign ificantly larger number of YSO candidates compared to the list of iRivera-Ingraham et al. 
(120111) . We suspect that the differences are mostly due to the additi on of X-ray and near- 


infrared criteria which increased the number of candidates in each class. iRivera-Ingraham et al. 


only list candidates with detections in all IRAC bands, which complicates a direct assessment 
of completeness. 


In Figure m we show the K and [3.6] band brightness distributions of the three identihed 
YSO classes in our sample. In the case of the [3.6] band, the three distributions show a sharp 
drop at a similar limit of about 15.5 mag. In contrast, for the K band, the distribution for 
Class III sources drops at about 15.0 mag, while the distributions for Class I and Class II 
sources drop near 18.0 mag. These histograms show that, due to the high and non-uniform 
extinction, the brightness distributions for young sources are likely incomplete above the 
sensitivity limits. However, our near-IR observations are deep enough to detect young sources 
across the entire region. For this reason, analysis of the luminosity functions discussed in 
Section 13.31 were performed using extinction-limited samples. 


3.2. Spatial Distribution of Young Stellar Sources 


In Figure [5] we illustrate the spatial distribution of YSO candidates in the W3 region, 
overlaid on an dust extincti on map construc ted with the optimized near-infrared excess 
method NICEST method of iLombardil (120091 ). The YSO sources appear to be very well 
constrained to the three well known areas of the complex (namely IC 1795, W3-Main and 
W3(OH)). 


Class I sources are mostly conhned to the areas of high extinction: following the pre¬ 
scription of Ihada et al.l (120131) . we found that over 80 percent of the Class I sources are 
located in high extinction regions {Ay > 7.0 mag), and the surface density of Class I sources 
is linearly correlated with Ay within 3 < Ay < 12 mag. This is in very good agreement 
with ot her cluster forming regions like Orion A (ILada et al.ll20131) and the Rosette Molecular 
Cloud ( Ybarra et al. 2013 ). 


Class II sources clearly permeate the entire complex, while Class III sources appear 
to have a slightly more centrally condensed distribution in the IC 1795 region but are also 
present in the western edge of W3-Main and surrounding the areas of low extinction to the 
East of W3(OH). 


Two aspects are important to notice. First, cluster regions are not separated clearly 
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Fig. 4.— Above: K band brightness distribution for young stellar source candidates identified 
in the W3 region. Below: [3.6] brightness distribution. 
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Fig. 5.— The plot shows the positions of infrared yonng stellar sources in the entire W3 
region. Red symbols indicate sources classified as Class I. Green symbols indicate Class II. 
Blue symbols indicate Class III sources. The background image and contours correspond 
to a visual extinction {Ay-, in units of magnitudes) map of the region constructed with the 
NICEST method. The map is Nyquist sampled over a 18" pixel grid, smoothed with a 
Gaussian hlter with a 36" FWHM. The contours indicate extinction at Ay =7.0, 10.0, 15.0, 
20.0 and 25.0 in units of magnitudes. 
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in the layout of the complex, which suggest at least some partial overlap along the line of 
sight perpendicular to the held. Second, the spatial distribution of young sources suggest 
that clusters are not single, centrally condensed systems. Instead, they appear to have 
sub-structure. 

In the following we a) make a consistent separation of the individual cluster populations 
and b) provide evidence of the presence of sub-structure in the cluster population by using 
a surface density analysis. We make the following two dehnitions: 

1. We dehne a Principal Cluster as a surface density structure that can be associated to 
a single Gaussian peak with a FWHM of at least 1.0 pc (see 1 13.2.11 below). 

2. We dehne a substructure as a signihcant surface density peak, smaller than a principal 
cluster. We identihed as signihcant, all those peaks with a surface density above 1.0 
stars per sq. arcmin (see 1 13.2.21 below). 


3.2.1. Identification of Principal Clusters 


In order to delimit individual clusters, we e mployed a Gaussian Mixture Model (GMM) 
analysis via the mclust p ackage on the R code flBanheld Sz Raftervlll993l: iFralev Sz Rafterv 


2002 


Rafterv et al.l 120121) . The package provides a model-based clustering analysis of any 


population, based on covariance parametrization and an Expectation-Minimization (EM) 
algorithm. It can also divide a population in an optimized number of clusters by means 
of the Bayesian Information Griterion (BIG). The BIG is the value of the maximized log- 
likelihood with a penalty on the number of model parameters, and allows comparison of 
models with diherent parameter sets as well as solutions with diherent numbers of clusters. 


We chose the Glass II sample as it is the most complete and it is distributed more or 
less uniformly across the complex. RA,DEG positions were fed into the code as a sample 
of bivariate observations. We chose to use circular Gaussians to minimize the number of 
parameters required for the estimation. Given the complicated layout of the clusters in W3, 
we opted to run the GMM analysis including noise modeling and rejection of outlier^. For 
our dataset, a model of hve clusters provided the best £t for the sample. This is shown in 


® In our case, we chose to use a first order estimation of noise in the spatial distribution of the sources 
on the plane based on nearest-neighbors. This way the mclust package applies the hierarchical clustering 
algorithm to a denoised dataset, and finally,reiects f rom the final solution, those sources that do not belong 


any cluster in the final solution (IRafterv et al. 


201211 . 
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Fig. 6.— Identification of principal clnsters in the W3 region nsing a Gaussian Model 
Mixture (GMM). This hve cluster model, applied to projected distribution of the population 
of Glass II YSO candidates in the region provides the optimal mixture based on a Bayesian 
Information Griterion (BIG). The circles are placed at the cluster centers determined by the 
GMM analysis and have radii equal to the Gaussian a. 
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Figure ini While this hnite mixture model shows less detail than -for instance- the nearest 
neighbor analysis, we see that it succeeds at identifying the three well known clusters in the 
region (W3-Main, W3(OH), IC 1795). It also adds two more principal clusters, one coinciding 
with NGC 896, the compact HII region located at the Western edge of the complex, and one 
more located at the Eastern edge. We labelled this group IC 1795-N. In Table |2] we list the 
result of the GMM analysis, with the centers and radii of the hve principal clusters identihed 
(columns 1 to 4). 


We emphasize that the purpose of this exercise is to dehne the principal clusters as 
consistent samples to study the history of star formation in the complex. The GMM analysis 
provides such consistent identihcation, whic h may differ with respect to other identihcation 
methods (e.g. 


Rivera-Ingraham et ahl 120111) . 


3.2.2. Cluster Substructure in W3. Nearest Neighbors. 


Estimating the surface density of young stellar sources is often done by applying a variant 


of th e nearest neighbor method (formalized for star cluster detection by ICasertano Sz Hut 


2010 


1985 ) . This method has been ap plied successfully in different studies (e.g. iGutermuth et ah 


Roman-Zhfiiga et ahll2008h . However, two very important issues regarding this method 


must be pointed out. First, when calculating distances to the nth neighbor in a distribution 
of points on a plane, there is a minimum value n = 6 to avoid biasing the detection toward 
small structures formed purely by chance in a 2-dimensional distribution of points. Second, 
when mapping the surface density dehned by such distances, the choice of an optimal smooth¬ 
ing kernel that would not dilute genuine substructures neither genera te spurious one s , is no t 
straightforward. These two issues were taken into consideration bv iGladwin et al.l (119991) . 
who developed two nearest neighbor mapping methods, gather and scatter, both optimized 
to correctly detect small groups (in our case, possible cluster substructures) in source posi¬ 
tion maps. The gather algorithm is more adequate for determining hierarchical structure, 
while the scatter method is more adequate for dehning the size of the smallest group in a 
distribution of points on a map. 


The gather algorithm was better suited for the task of detecting and separating sub¬ 
structures. We applied gather to the distribution of all Glass II YSO candidates in the W3 
region, which is the most complete of our samples. Densities were calculated with a square 
kernel of 36'' over a Nyquist sampled RA-DEG grid. In Figure [7] we show the resultant map, 
in the form of surface density contours overlaid on an infrared emission image that combines 
data from the Spitzer MIPS 24/im, 70/im from Herschel PAGS and 250 pm from Herschel 
SPIRE. If we compare the density map to the spatial distribution of sources in Figure [5] 

























17 



Fig. 7.— RGB image constructed by combining map from Spitzer MIPS 24/im (blue), Her- 
schel PACS 70 pm (green) and Herschel SPIRE 250pm (red). Contours indicate surface 
density of Class II young stellar sources the W3 region, estimated with the Gather method. 
Contour levels are set as 0.5,1.0,2.0,3.0,4.0,5.0 and 10.0 stars per square arcmin,. The posi¬ 
tions of Class I sources are indicated with red circle symbols. The positions of O and B type 
stars are indicated with magenta and cyan color symbols, respectively. Areas in yellow/green 
are saturated pixels in the infrared images. The labels correspond to previously identihed 
clusters and sub-clusters in the W3 region ; labels starting with letters “RI” indicate those 
identihed bv iRivera-Ingraham et al.l (120111) . 
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we see how gather highlights the three main regions (IC 1795, W3-Main and W3(OH)) 
as the principal clusters, but also shows a number of peaks that clearly dehne signihcant 
sub-structure in the spatial distribution of young sources. 

The positions and equivalent sizes of sub-structure peaks in the gather map are listed in 
Table [H Some of the sub-structure peaks in the map coincide in position with one or several 
Class I sources. Those peaks are likely related to the most recent episodes of formation in 
the complex. Notice how none of the peaks near the center, corresponding to IC 1795, are 
associated with Class I sources, possibly indicating a n absence of current formation. We also 
show the p ositions of O ari d B type stars listed by lOgura fc Ishidal (jl976|), iNavarete et ah 
( 2011 1 and Bik et ah ( 2012 1. 


The largest concentration of massive stars (6 0-type stars, 7 B-type stars) is located at 
W3-Main, along with a noticeable concentration of protostars, associated with at least four 
YSO surface density sub-structures. In IC 1795 there are tw o 0-type stars and one B-type 
star identihed. The central O star in IC 1795 (No. 89 in the lOeura fc Ishidal (119761) list) is 
associated with the highest surface density peak. In the W3(OH) region, two sub-structures 
dehne the positions of groups B and C in W3(OH), and both are associated with B type 
stars. It is interesting that the main cluster (group A) in W3(OH) is not associated with 
any identihed massive source. 













No 

1 

2 

3 

4 

5 

6 

7 
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9 

10 

11 

12 
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14 

15 

16 

17 

18 

19 

20 

21 

22 

23 

24 

25 

26 

27 

28 

29 


- 19 - 


Table 1. Sub-structures in W3 


Center Coords. Class P Notes 
J2000 


02:27:33.29 

61:56:24.2 



02:27:21.43 

61:54:57.4 


W3(OH)-C 

02:27:16.37 

61:54:16.0 


W3(OH)-B 

02:27:12.43 

61:53:08.8 

1 


02:27:03.46 

61:52:33.3 

3 

W3(OH)-A 

02:27:18.20 

62:00:04.1 

1 

R106 

02:27:08.58 

61:56:30.6 

1 


02:26:58.74 

61:54:53.8 

4 


02:26:49.52 

61:57:46.0 

1 


02:26:48.17 

62:03:28.3 



02:26:40.22 

62:00:18.4 


IC 1795 center 

02:26:41.14 

62:04:33.6 


R107 

02:26:38.87 

62:02:01.3 



02:26:38.87 

62:00:42.2 



02:26:35.50 

61:59:15.1 



02:26:33.25 

61:56:25.0 



02:26:20.01 

62:05:01.3 

1 


02:26:16.96 

61:59:48.7 



02:26:19.22 

61:56:42.8 


RlOl 

02:26:03.36 

62:07:13.7 



02:25:59.28 

61:58:37.3 



02:25:57.97 

62:08:50.6 

1 


02:25:52.42 

62:04:21.4 

2 

RlOO 

02:25:51.10 

62:00:20.0 



02:25:47.76 

62:10:35.2 


R104 

02:25:42.74 

62:08:10.7 



02:25:43.36 

62:06:18.0 


W3-Main center 

02:25:32.91 

62:06:49.4 



02:25:32.98 

62:04:32.9 

2 
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In Figu re [7| we also placed labe ls (“RI”) next to regions coinciding with “sub-clusters” 
identified by iRivera-Ingraham et al.l . It is worth noticing how the gather map highlights the 
groups A, B, C that were identihed by eye in Figure [21 


3.2.3. YSO Ratio Maps 


It is possible to trace the evolution of a star-forming region through the construction of 
maps of YSO number ratios. Figure [8] shows maps for the Class II to Class III number ratio, 
Rn:iib and t he Class I/O t o Clas s II number ratio, Ri-ji, both constructed using the method 
described in [Ybarra et ah! (j2013|). 


For the Rilih map we limited the source candidates to a de-reddened brightness H=13.5 
mag. This brightness cut avoids biasing the number counts fro m incompleteness due to 
patchy extinction. We used the Bayesian estimator described in [Ybarra et ahl (120131) to 
calculate the ratios within circular projected regions with radius 3'. The map only shows 
ratios for regions with sufficient number counts such that iVi -|- W -|- 1 < 0.25(iVi -|- l)(iV 2 — 
1), where Ni is the number count of the numerator, and W is the number count of the 
denominator. 


We found that the whole complex has Rilih >1.0 which suggests star formation was 
ubiquitous throughout the complex ~ 2-3 Myr ago. We can also trace the most recent star 
formation by using the Rlh ratio (Fig. [HI right). Regions with signihcant are expected 
to have been actively forming stars within the last Myr. We found that the regions of high 
Ri:!! appear to form an oblong ring around IC 1795. This mirrors the distribution of Class 
I/O sources which also appear to surround IC 1795 (see Fig. [7]). Interestingly, the high Rilih 
in IC 1795 suggests that stars were forming there up to 1 Myr ago. Therefore the absence of 
Class I/O sources in IC 1795 suggests that star formation was recently and possibly abruptly 
terminated in that region. 

It is important to notice that values for the Rilih ratio in the map can be high (typi¬ 
cally larger than 5), particularly because the Class III sources are more spatially scattered. 
However, the Rimi is not the same as the circumstellar disk fraction in the cluster, because 
we are only measuring against weak disk sources and not the whole population (see also 
section [H751 below). 
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Fig. 8.— Top-. Class II to Class III ratio map. The black contours indicate = 1.0,3.0. 
Red circles indicate the positions and radii of principal clusters as in Fig. [6l Bottom-. Class 
I/O to Class II ratio map. 
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3.3. K band Luminosity Function of Clusters in the W3 Complex 


We obtained the K band luminosity function (KLF) for each of the hve principal clusters 
identihed with our GMM analysis using the following procedure: 

Each cluster sample is composed by all K-band sources detected within a “3-sigma” 
radius dehned as i? 3 o- = (3.0/\/81n2)rci, where is dehned from the GMM analysis as 
the half-width at half-maximum of the best ht 2-D Gaussian. Dehning samples within 
minimizes effects of cluster overlapping, which could r nislead in the estima tion of age and 
age spreads in cluster forming regions (e.g. NGG 1980. lAlves fc Bouvll2012l) . 


We obtained an extinction limited sample, dehned as follows: 1) we pro jected a Ay =20 
mag extinction vector from a 1 Myr isochrone flD’Antona fc Mazzitelhlll997|) until it reached 
the sensitivity limit of the observations. The Ay-limited sample is dehned by all sources 
enclosed within the area dehned by the extinction vector, the zero-age main sequence and the 
reddened isochrone (see Figure IHl left). 2) we separated the possible foreground population 
from sources with colors coinciding with the locus of the dwarf and giant sequences in the 
J — H vs H — K color-color diagram.This dehnes a foreground contribution. In Figure 
[9] (right panel) we show a color-magnitude diagram for W3-Main showing the Ay-limited 
sample selection. 


We estimated the fraction of sources with circumstellar disks within R^cr, for each prin¬ 
cipal cluster, as the total number of Glass I plus Glass II sources in the extinction limited 
sample, down to a de-reddened brightness H=13.5 mag. (see also section [3.2.31 above). The 
highest fractions of circumstellar disks, close to 30 percent, are in W3-Main, IG 1795 and 
W3(OH). Glusters NGG 896 and IG1795-N have smaller fractions, closer to 15 percent. 

The total contribution from the held was estimated from a nearby control held with min¬ 
imum extinction, observed in the same conditions and to the same depth as the cloud helds 
(see Appendix [^. We made a K-band histogram from the control held catalog and scaled 
it to the same area of the cluster sample. Then we subtracted the foreground contribution 
from the control held KLF, and applied artihcial reddening to the remaining distribution 
down to the same level as the cloud sample. Artihcial reddening is applied by constructing 
a normalized Ay extinction distribution from the background stars in each cluster region 
using the NIGER algorithm, and used it as a frequency distribution to convolve the control 
held K-band histogram, producing a hnal control held K band brightness distribution. The 
NIGER algorithm provides estimates of extinction with average errors aAy <0.1 mag be¬ 
low Ay =30 mag. These errors are way below the size of the K-band brightness bins used 
to construct the KLF. The hnal KLF of the cluster sample is obtained by subtracting the 
hnal Ay-convolved (reddened) control held distribution from the one constructed from the 
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K-band limited sample. 


The advantages of using a K-band limited sample to construct the KLF are: a) that 
all samples are cut at common maximum extinction; b) that the contribution of background 
extra-galactic sources is reduced to a minimum; c) that patchy extinction is less prone to 
bias the brightness distribution, specially at the most embedded regions -e.g. centers of 
W3-Main and W3(OH). Also, by convolving the raw control held KLF by the “on” held 
extinction distribution, we are able to obtain a fair estimation of the contribution from held 
stars on diherent lines of sight flMuench et al.l 120031) . 


3.3.1. Age Inference from Cluster Model Fits 


The KLF for each principal cluster sample was then compared to artificial K-band 
brightness distributions of embed ded cluster po p ulatio ns of varying ages, constructed with 
the model interpolation code of iMuench et al.l (120001) . The method is relatively easy to 
apply: using a set of pre-main sequence models, the code uses a Monte Carlo generator 
to draw an artihcial population from an assumed initial Mass Function (IMF) with a total 
number of members input by the user, between a minimum and a maximum ^e, dehned 
by model isochrones, properly shi fted by the distance modulus of the complejo. For this 
study, we chose the isochrone set of iD’Antona fc Mazzitellil fjl997|) . with a deuterium fraction 
[D/H] = 2 X 10“®. The user has the option to choose the form of the IMF and also has 
the option to hx or to sample values for a set of parameters that includes an extinction 
distribution, a binary fraction, and the infrared excess fraction as a function of color (H — K). 
The peak and the slopes of the KLF are particularly sensitive to age and age spread, which 
are the two parameters that were left to vary. We decided to l eave t he binary fraction as 
a hxed, value of 20%, based on estimations by IPadgett et al.l (119971) suggesting that the 
binary fraction in the Orion B clusters may be about 15% in separatio ns up to 1000 AU 
-the separations we are able to resolve in our near-IR images. Also, iMuenchl (120021) in 
a similar experiment showed that, particularly, the mean age and spread estimates were 
not signihcantly different when changing the binary fraction within 40 percent. A smaller 
fraction is probably too optimistic, but a fraction that is signihcantly large (30-40%) tends 
to add dispersion in the histogram and “dilutes” the location of the KLF peak. For the IMF, 
we decided to use the the parameters of the 3-part broken power law obtained for the IC 348 


® The 0.107 kpc distance error quoted by Hachisuka (2006) for their maser parallax distances, translates 
into an error of about 0.2 mag in the brightness estimates. This is smaller than the size of our KLF bins 
(0.5 mag.) In terms of mass, this translates into an uncertainty of about 0.05 in a log-mass scale. 
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cluster by iMuench et al.l fl2003fl l^u. Each of the hxed parameters will add its own error to the 
estimate of the best value for the running parameters, and there are also errors associated 
with the use of one or other isochrone model set. A full discussion of the method is out of 
the scope o f this p aper, however, the feasibility of this method has been properly discussed 
by Muench ( 2002 1 (chapter 5). 


We made a grid of cluster models that varied the minimum age of the population between 
0.5 and 10 Myr, in steps of 0.5 Myr, and having star formation age spreads of 1.0 to 5.0 Myr. 
For each minimum age and age spread value, we simulated 500 artihcial KLFs and compared 
each of them with our observed KLF using a simple reduced scheme. Averaging over the 
simulations, we are able to determine which age or group of ages hts the observed KLF better 
(see Figure in]). Our estimate of the age and age spread is dehned by those models for which 
the conhdence limit is above 68% (2-sigma). In Figure ITU] we show two examples of the age 
estimation scheme: we plot the mean age used in the simulation versus the width of the age 
spread. The reduced value contours indicate the grid values for which the age models ht 
better (25 percentile). In these examples, the W3-Main cluster shows a very clear ht for a 
population with an age spread of less than 2.5 Myr, and a mean age of around 2.0 Myr. The 
IC 1795, instead, is less well constrained: we see two local x^ minima in the contour map. 
The hrst one agrees with a mean cluster age of less than 2.5 Myr old and an age spread of 
less than 3 Myr. The second minimum suggests a mean age between 2.5 and 3.5 Myr and an 
age spread between 3 and 5 Myr. The sec ond option agree s much better with spectroscopic 
studies of IC 1795, particularly the one bv lOev et al.l (120051) . To be conservative, we consider 
a mean age between 2.0 and 4.0 Myr, and an age spread in the whole range of 1.0 to 5.0 Myr. 
The main reason for the discrepancy in IC 1795 is the large number of very young (Class 
II) sources still present in the IC 1795 region, which is indicative of a mixed age population. 
The analysis suggests a very similar age and spread for W3(OH), which is surprising given 
that it hosts more embedded populations. The analysis for clusters IC 1795-N and NGC 896 
yield mean ages from 2.5 to 4.5 and from 3.5 to 5.0 Myr with possible spreads from 2.0 to 5.0 
and from 1.0 to 5.0 Myr, respectively. Unfortunately, the upper limits in the age spread are 
not well constrained for most clusters, however, age spreads much larger than 5 Myr sound 
little plausible for embedded star cluster populations. In table |2] we list the main results of 


^^Notice that we cannot sample the whole mass spectrum at the distance of W3, specially given the high 
extinction. However, our Ty-limited samples are complete down to K=16.25 mag which at the distance of 
W3 (2.05 kpc) and given the estimated ages of the clusters, could range between 0.05 and 0.3 M 0 , with an 
average around or below 0.1 M©. This is below the expected peak of the IMF at 0.2 M© and that is why 
we con sider necessary to us e the three parts of the broken power-law parametrization. The interpolation 
code of Muench et al. ( 2000[) takes into consideration the minimum and maximum brightness of the cluster 
observations and samples the IMF within the proper ranges. 





















25 


10 - 


15 - 


• ■' a?ofr2 ^ O no ° 



20 


- 0.5 0.0 0.5 1.0 1.5 2.0 2.5 

J-K [mag] 



Fig. 9.— Left: Extinction (Ay) limited selected sample for W3-Main using K vs H — K 
color-magnitude diagram. The solid thick line is the z ero-age main sequence. The tw o thick, 
dotted lines are a 3 Myr isochrone from the models of iD’Antona fc Mazzitellil (119971) . before 
and after adding an extinction of Ay = 20 mag. This is also indicated by the extinction 
vector that runs from the lowest mass point (0.3 Mq). The dot-dashed line indicates the 
sensitivity limit of our Calar Alto data set. Solid symbols indicate sources selected out as 
possible foreground sources. The open circle symbols indicate sources selected as cluster 
members in the hnal sample. Right: An illustration of KLF htting, for IC 1795. The solid, 
black histogram is the observed KLF, after subtraction of the held component; the red, green 
a nd blue lines ar e artihcial KLFs modeled after the properties of the cluster, using the code 
of iMuench et al.l. with reduced probabilities of 0.95, 0.68 and 0.22, respectively. 
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Fig. 10.— Age estimates for all principal clusters, illustrated with a contour map of the 
normalized values. We plot the mean age input for the artificial KLF models vs. the age 
spread used in each simulation set. Areas at and below purple color levels indicate the most 
likely values. The white contour lines trace the 68 and 95 percent conhdence limits. The 
light line surrounding the contour areas indicates the limits of the model grid. 
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the artificial KLF modeling analysis for the hve principal clusters. The age and age spread 
for each of the principal cluster is listed as the center and width of the top contour in the 
age spread vs age maps. 


3.4. The Context of the W3 Molecular Cloud 


The 1.1 mm Bolocam maps and the CO isotope radio emission maps of iBieging fc Peters 


(120111) show that the IC 1795 is almost free of molecular gas emission, while most of the 
remaining cloud is concentrated on two large clump systems at W3-Main and W3(OH). The 
latter has a total mass of approximately 1.9 x 10^ Mq, while the former has a total mass 
of approximately 2.5 x 10^ M ^t^, divided in two main clumps named W3-East and W3-West 
(iRivera-Ingraham et ahl 120131) . In such layout, IC 1795 is seen surrounded by an almost 
symmetric shell of star forming molecular material. 


In Figures ITT] and \n\ we show the locations of Class I and Class II source candidates 
in the context of the ^^CO(l-O) hrst moment (average radial velocity) map of the cloud, as 
well as a map of 1.1 mm emission over infrared images at both the Eastern (W3(OH)) and 
Western (W3-Main) sides of the Cloud. 


The W3-Main map shows a very notorious blue-shifted peak at the center of the principal 
cluster, that coincides in velocity with a hlamentary structure that extends south toward 
NGC 896 and north toward IC 1795-N. The blue-purple colored peak corresponds to a 
velocity range centered around -42.5 km/s, and it ha s been suggested to ind icate central 
massive outflow activity at the center of the cluster flBieging fc PetersI l201ll) . The main 
^^CO(2-l) emission structure associated with W3-Main with velocity range centered around 
-40 km/s (green colored in the image), shows a clumpy structure suggestive of other local 
velocity gradients. A few small receding clumps associated with velocities in excess of -36 
km/s (red colored in the image) possibly represents gas still being expelled from the central 
cluster, and conhrms that gas removal in IC 1795 was very recent. 


The W3(OH) map shows that the cluster formation activity at this side of the complex 
is associated with molecular emission with signihcant radial velocity gradients. This eastern 
side of the cloud has a hlamentary morphology, with almost uniform radial velocity. There 
is a 2 km/s difference between the east and the west edges of the hlament, suggestive of 
a “rolling” movement of gas in the north-south direction, but on the eastern side the peak 
velocity alternates several times from 46 to 48 km/s, indicating local gradients possibly 
linked to gas compression and expulsion at the W3(OH) groups. This is particularly evident 
at group W3(OH) A. 
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Table 1—Continued 


No. 

Center Coords. 

J2000 

Class P 

Notes 

30 

02:25:30.88 

62:00:05.7 

2 

NGC 896 center 

31 

02:25:30.41 

61:57:11.6 

3 

R103 

32 

02:25:21.90 

62:07:06.8 



33 

02:25:10.23 

62:03:11.0 


R105 

34 

02:25:13.21 

61:59:21.6 



35 

02:25:00.93 

62:09:08.7 

3 

R102 

36 

02:25:03.58 

62:06:52.4 

1 


37 

02:25:00.26 

62:05:37.0 

3 



’^Number of Class I candidates within gather surface den¬ 
sity peak contour 


Table 2. Main Parameters and Age Estimates for Principal Clusters in W3 


Cluster 

Center Coords. 

Tcl 

Disk Fraction’’ 

Mean Age'^ 

Age Spread'^ 


J2000 

[pc]®" 

[%] 

[Myr] 

[Myr] 

W3-Main 

36.368671 

62.097485 

1.313 

27±4 

2.0 [1.5,2.5] 

1.5 [1.0,2.5] 

IC 1795 

36.652903 

61.990908 

1.244 

29±5 

3.0 [2.0,4.0] 

4.0 [1.0,5.0] 

W3(OH) 

36.798850 

61.895349 

1.644 

27±3 

3.0 [2.5,3.5] 

3.5 [3.0,5.0] 

NGC 896 

36.388061 

61.996607 

1.428 

14±3 

3.5 [2.5,4.5] 

4.0 [2.0,5.0] 

IC 1795-N 

36.569302 

62.103874 

1.618 

12±2 

4.0 [3.5,5.0] 

2.0 [1.0,5.0] 


'^Assuming a distance d = 2.045 kpc to the W3 Complex. 

’’Estimated as the fraction of Class I+II sources within i? 3 o-, in the extinction-limited 
sample, down to Hdered=^3-5 mag. 

^The number pairs inside the brackets indicate the range of model ages within the 2- 
sigma (68%) conhdence level. The number listed before the bracket is the central value, 
which we list as the best estimate for the age/age spread of the cluster 
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Fig. 11.— Association of gas and stars in the W3-Main region. Left: The RGB 
image in the backgronnd is constrncted from Spitzer-IRAC [3.6],[5.0] and [8.0]pm mo¬ 
saics, and the contonrs indicate 1.1 mm emission in Jy/beam from the Bolocam map at 
0.12,0.24,0.48,0-96,1.2,2.4,4 -8 and 9.6. Right : K hr st moment ^^CO(2-l) map of the same 
region, from the dataset of iBieging fc PetersI (120111) . The color table indicates radial veloc¬ 
ity in km/s. The red and orange dot symbols indicate the positions of Class I and Class II 
sonrces identihed from onr main catalog, respectively. 
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Fig. 12.— Association of gas and stars in the W3(OH) region. Symbols and colors are 
equivalent as those in Fig. [11] 
















31 


In Figure [T^ we show two position-velocity (p-v) plots corresponding to two cuts along 
the main structures of W3(0H) (cut LI) and W3-Main (cut L 2) in the 1.1 mm r nap. We 
constructed these p-v plots from the ^^CO(3-2) emission map of iBieging fc PetersI (120Ilf) . 


The LI plot shows a well dehned bipolar structure at the location of cluster A with 
a peak to peak velocity difference of 15-20 km/s within less than 0.5 pc. The ^^CO(3-2) 
map does not resolve the structure well enough to claim the presence of individual outflows, 
but the plot shows that the W3(OH) elongated clump may be associated with a structure 
coherent within (roughly) -50 and -45 km/s. The plot shows an opening flow of material 
associated with cluster A (offset~ 170"), and a tighter blue-shifted flow associated with 
cluster B (offset~ 200"). A third component may be present near offset~ 300", possibly 
associated with cluster C. 


The L2 plot shows a very strong outflow feature associated with the main molecular gas 
clump in W3-Main, with a peak to peak velocity difference of 25-30 km/s within less than 
0.7 pc. The red-shifted component is associated with the left lobe of the main molecular 
clump and vice versa. The two peaks of the red-shifted structure suggest a parabolic opening 
indicative of a widening outflow. The physical scales and velocity differences of these struc- 
tures are very similar to those observe d in other mas sive star for ming regions like DR-21 


(Schneider et ah 

2010 

), G240.31+0.07 ( 

Oiu et ah 

2009 

) or W49A ( 

Peng et al. 

2010) 


We estimated the integrated intensity of ^^CO(3-2) in two square areas of 120" x 120" 
at the center of the double peaked clump in W3-Main, which correspond to the blue and 
red peaks of the outflow feature. We deh ned these square areas from a hrst moment (mean 
velocity) map. Using the prescription of iGinsburg et al.l (1201 ll) . we estimated the column 
density, N{H 2 ) and the total mass of outhow gas, Mout enclosed in each th e two area s (17.6 
and 40.7 Mq for the red and the blue peak, respectively). Then, following iQiu et al.l (1200911 
we estimated the dynamical age, t^yn of each component (6.5 and OxlO'^ yr, for the red and 
blue peaks, respectively) and the mass loss rate, Mout = Mout/tdyn- We estimated that the 
total mass loss rate is close to 7.2 x 10 “^Mq yr“^. This is close to estimates in other massive 
star forming regions, which have mass loss rates of ~ 10 “^Mq yr“^. A mass loss rate like 
this, if constant, could remove 10^ Mq of gas in about 1.5 Myr. 


4. Discussion 

Our KLF analysis suggests that episodes of cluster formation in W3 may have started 
over 4 Myr ago, without any well dehned pause. Actually, our data also suggests that star 
formation has been ubiquitous across the whole complex within the last 1.0-2.0 Myr: 
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Fig. 13.— Top: Bolocam 1.1 mm emission map (grayscale, overlaid with same set of contours 
as in Fig. m- Arrows indicate length and direction of two cuts used to create position- 
ve locity maps. Midd l e and Bottom: Position-velocity plots, from the CO(3-2) emission maps 
of iBieging fc PetersI (120111 ). The color table indicates line temperature in K. 
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1. All of the principal clusters exhibit a signihcant fraction of sources with detectable 
circumstellar emission. This is evidence in favor of recent formation across the whole 
complex. Class I sources are distributed preferentially in high extinction regions, in¬ 
dicating the most recent episodes of formation. However, Class III sources are present 
also near these highly embedded regions: this could be indicative of age spreads com¬ 
parable to the T Tauri timescale (and possibly longer, if we consider young sources that 
escape classihcation via our color and X-ray emission scheme) across the whole com¬ 
plex. However, this could also be the result of rapid dissipation of disks, particularly 
for intermediate to high mass sources. 

2. The age spread of IC 1795 cluster may be at least 4 Myr. This cluster also has a slightly 
larger density of Class HI sources (e.g. 7 vs 5 percent of all YSO candidates down to a 
de-reddened brightness of H=13.5 mag, estimated within in IC 1795 and W3-Main, 
respectively). This is in agreement with previous studies that suggest this is the more 
evolved population in the complex. Also, the lack of nebulosity or strong CO emission 
is indicative of gas dispersal, which agrees with a more evolved cluster. However, the 
mean age that we estimate for this cluster is not very different from that of the youngest 
clusters in the complex, W3-Main and W3(OH). Also, the Rn-.iii ratio in IC 1795 is 
comparable to W3-Main and W3(OH) and so is the total fraction of circumstellar disk 
sources in this cluster (about 30 percent within r3o-[IC1795]). Still, no Class I sources 
are longer present in IC 1795. These results are not contradictory. They all suggest 
that the forming gas dispersal timescale in IC1795 was rapid, comparable or faster 
than the T Tauri timescale. 

3. The KLF analysis suggest that W3(OH) has an age and age spread not very different 
from IC 1795. Again, this is indicative of more than one episode of formation. The 
main groups A, B and C in W3(OH) are still embedded in bright nebulosity and the 
principal cluster has a large circumstellar disk fraction, comparable to IC 1795 and 
W3-Main. There is also a signihcant number of Class I sources in W3(OH). However, 
there may be slightly older groups toward the eastern edge of the complex, judging 
from the presence of Class H and Class HI sources to the edge of the molecular clump 
and no nebulosity. The Rhjh maps show an abrupt termination at the eastern side 
of the W3(OH) principal cluster which reinforces this idea. The structure in W3(OH) 
(specially if there is an older component) weakens the argument of a Gaussian cluster 
used for the GMM analysys. We need deeper photometric data in order to reconcile 
this problem. 

4. The KLF modeling analysis shows that IC 1795-N and NGC 896 may host the oldest 
clusters in the complex with a mean age between 3.5 and 4 Myr. These two aggregations 
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could have been formed in an external layer of the original molecular cloud, flanking 
IC 1795 from north and south. 

5. Our analysis suggests that W3-Main hosts the most recent cluster formation episode, 
starting just over 2.0 Myr ago and having the shortest spread, with less than 2.5 Myr. 


In Figure [m we present a cartoon illustrating one possible scenario of the progression of 
cluster formation in W3 based on our KLF analysis. We suppose that the original molecular 
cloud might have been an area somewhat larger than the current span of the cluster complex, 
over 4 Myr ago. As shown bv iThronson et ahl fjl985l) . the W3 complex formed in a ridge of 
swept up molecular gas in the western age of the giant HII region W4. The expansion of W4 
may have triggered the collapse of the W3 molecular ridge. By then, the formation of IC 
1795-N and NGC 896 might have started, while IC 1795 and an area of unknown extension 
near the eastern edge of W3(OH) (we indicate this region symbolically in the cartoon) started 
to collapse and form stars. Just over 3 Myr ago, the process of cluster formation continued 
with IC 1795 and the hrst episode in W3-Main and W3(OH), while clumps corresponding 
to the present W3(OH) cluster group would be collapsing. About 2-3 Myr ago, gas dispersal 
could have started in IC 1795, IC 1795-N, NGC 896 and the eastern edge of W3(OH). The 
most recent episodes in W3-Main and W3(OH) would have started. By then, gas would be 
dispersed away from the center of the complex and start concentrating at the W3(OH) and 
W3-Main regions. This would resemble the current layout of the complex, with most of the 
remnant molecular gas being at W3(OH) and W3-Main. In this scenario, cluster formation 
progressed hrst from north and south towards the center of the complex and then from the 
center towards east and west.It is important to notice that this scenario is based on our 
estimates of the mean ages of the clusters, but could only be supported by a better constrain 
of the age spreads of the clusters. 

One important point to discuss here is whether or not induced formation plays a domi¬ 
nant role on the formation of cl usters in the W3 complex, as suggested by previous studies. 


For instance, iRuch et al.l (120071) support the picture in which IC 1795 trigge r ed th e forma¬ 


tion of W3 Main and W3(OH) episodes. The study bv iRivera-Ingraham et al.l (120111 ) suggest 


that triggering processes were important in the formation of the clusters hanking IC 1795, 
but also that triggering processes work at local (sub-parsec) scales, with high mass stars 
acting to conhne and compress material, enhancing the efficiency for the formation of new 
high mass stars by making convergent hows. This, in turn, would help to the progression 
of formation from the oute r toward in the inner regions of the complex. The study by 


Feigelson fc TownslevI (120081) . pointed to signihcant gradient of ages between the central OB 


sources of the cluster and a widely distributed population of PMS stars. 
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What does our study tell about the global star forming history of the complex? On 
one hand, if gas removal in the central cluster IC 1795 contributed to collect gas in a 
shell like structure, it may have enhanced the efficiency of formation in W3-Main and 
W3(OH). This would be in agreement with the flow convergence scenario discussed by 


Rivera-Ingraham et al.l . According to that picture, IC 1795 could have been the major 


source of induction for star formation in the molecular shell. 

Our KLF model analysis indicates large age spreads across the complex, with a progres¬ 
sion that lead to cluster formation across the entire complex in a period of about 3 million 
years. This is in good agreement with our analysis of the YSO class surface density ratio 
across the complex. In that case, the progression is possibly related to primordial density 
structure of the cloud. This is in agreement with the scenar i o proposed for other complexes 


like t he Rosette Molecular Cloud flRoman- Zuniga et al.ll2008l: IWang et al.ll2009l: lYbarra et ah 


20131 ). where it has been shown that progression of cluster formation could rather follow the 
density structure of the cloud than a process of induction in a sequence from the older to 
the younger episodes. 

Our YSO surface density analysis provides evidence of sub-structure in the principal 
cluster populations. Small aggregations of stars around the main cluster events (very ev¬ 
ident in W3(OII), and groups like R102 and R103), revealed as significant peaks in the 
surface density of YSOs, in many cases coincident with the locations of Class O/I candidates. 
This is suggestive of a complex scenario, in which stellar collections of different sizes may 
form with small age differences by the compression of gas in various locations within large 
clumps. Strong evidence o f sub-c lustering in W3-Main was also provided by the very deep 
near-IR KLF by iBik et al.l (120141) : they pointed out that an age gradient was evident from 
edge to center in W3-Main, judging from a higher fraction of circumstellar disks near the 
center. The large number of disk bearing sources in IC 1795 could be indicative of a similar 
process, where cluster formation occurs in sub-structures. If these sub-structures do not 
form simultaneously, they certainly may lead to a large age spread, as observed. 

The projected separations between surface density peaks in the gather map range is 
1.4±0.4 pc. These separations are comparable to the estimated radii of the principal clusters 
identihed with our GMM model. This scale is shorter than the Jeans length (2.0 pc) for a 
mean density of lO^cm^ and Tk = 10 K, typical of the diffuse gas traced by ^^CO and may 
be tracing the primordial fragmentation of the cloud. One possibility is that density peaks 
in the gather map are sub-clusters that eventually merge into larger groups, but that will 
depend on the dynamical stability of the aggregations an d the gravitational potential of the 


corre sponding principal cluster. Numerical simulations flBonnell et al.l 1201 ll: lAllison et ah 


20101) suggest, that clusters are initially highly structured, but rapidly evolve by merging 






















36 


c;-fl (tyr JQ* 


jr--\ ' 


CT ^ 1 

/ 

iliihn 

V • - " V 1 

V \ 

It,''’' 



Msri- 'JihJ 

t-fl? uyr ^11 


• n 

‘ 'v 

Jl. Bi- ’ 

K~i M4i f 

* Jilf kill? 1 



Fig. 14.— A cartoon illustrating our proposed scheme for the history of star formation 
in the W3 Complex. The blue regions indicate the extension of the molecular gas cloud. 
Dotted, empty regions indicate regions of collapse towards cluster formation. Red shaded 
areas indicate active cluster formation regions. Green shaded areas indicate cluster regions 
where gas dissipation has started and yellow shaded areas indicate slightly older clusters 
where gas dispersal has removed a signihcant amount of the parental gas and the cluster. 
The rectangular, easternmost area, represents a region of earlier formation with unknown 
extension, that might contribute today to a larger age spread observed in the W3(OH) 
cluster complex. Black arrows indicate the suggested progression of cluster formation in the 
complex. 
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into large entities. On the other hand, small gronps of stars may disperse rapidly or even 
remain stable depending crucially on the nu mber of members and the star forming efficiency 
flLada fc Ladall2003l: lAdams fc Mversll200ll) . It is also expected that most small, embedded 
aggregations end up being disrupted early, as proposed by the “infant mortality” and “cruel 
cradle” models. 


The formation of small groups that disperse rapidly and merge into a larger central 
structure could provide a simple explanation for the presence of larger age spreads in a lower 
surface density population in regions like IC 1795-N or NGC 896, and would be in agreement 
with age spreads and sub-structure in the younger principal clusters. 

The strong outflows in the molecular gas shell at W3-Main and W3(OH) conhrm a very 
rapid removal of gas in the complex. This removal may occur in a period comparable to the 
T Tauri timescale, as evidenced by Class II source candidates in IC 1795, where gas is mostly 
dispersed. However, it is important to stress that the mass loss rate cannot be constant, 
and will probably decline within a few 10^ Myr, if it has not already done so. The estimate 
we provided, however, may not be far from the average mass loss rate during the main gas 
dissipation period. 


5. Conclusions 

We obtained deep, high resolution ground-based near-IR images of the Eastern edge of 
the W3 region in the Perseus Arm. These data were combined with archive imaging data 
from the Spitzer and Chandra space telescopes to obtain a new catalog of young sources for 
the W3 Complex. We used our resultant catalog to: 

1. Identify and classify young (Class I, II, III) embedded sources in the region. 

2. Determine the locations and extension of the principal clusters in the complex. 

3. Determine the spatial distribution of young sources across the complex. This distribu¬ 
tion shows signihcant substructure. 

4. Construct extinction limited samples that allowed us to obtain clean K-band luminosity 
functions of the principal clusters. We compared the resultant functions with those 
from artihcial clusters of ages created from pre-main sequence models. This analysis 
permitted to estimate mean ages and possible age spreads for the principal clusters of 
the W3 Complex. 
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We found that IC 1795, the central cluster of the complex, still has a large fraction of 
Class II sources, despite not being embedded anymore, like clusters W3-Main and W3(OH). 
The Rn-.iii ratio is large across IC 1795, just as it is in the flanking embedded clusters. 
The ^^CO map shows some receding molecular gas clumps coincident with the center of IC 
1795. This, along with the absence of Class I sources suggest that parental gas was removed 
quickly from IC 1795, in a period shorter than the T Tauri timescale. The high mass loss 
rate estimate at the center of W3-Main supports this strongly. 

Our Gaussian Mixture Model analysis clearly identihed the well known clusters IC 1795, 
W3-Main and W3(OH). In addition, two principal clusters could be identified, one associated 
with the compact HII region NGC 896 and one more located north of IC 1795. 

We estimated the ages and age spreads of the principal clusters by comparing their K- 
band luminosity functions with those of artificial young clusters. The differences between the 
ages of the clusters are not signihcantly large. We also were able to constrain the age spread 
of the W3 cluster, but we could not constrain well the age spread of the other four clusters, 
so there might be overlaps in the formation periods of all clusters in the complex. This, is 
conhrmed by a large number of Class II and Class III sources across the whole complex. 
This suggest a progression of cluster formation scenario in which formation proceeded as a 
result of primordial density structure of the cloud. 

Using our list of Class II candidates, we constructed a surface density map that shows 
a signihcant amount of substructure in the principal clusters of the complex. The large age 
spread in the clusters may indicate that formation proceeds in episodes. If sub-structures 
observed do not form simultaneously, they certainly may lead to a large age spread, like we 
observe. 


We propose a scenario in which star cluster formation progressed over 4 Myr ago from 
north and south (clusters NGC 896 and IC 1795-N) toward IC 1795. Later (about 2 to 3 Myr 
ago), the progression moved from center (IC 1795) towards East and West (W3(OH) and 
W3-Main). We also suggest that an undefined region East of W3(OH) could form groups of 
stars before the formation of the embedded multiple embedded groups we see today. 


The formation of W3(OH) and W3 Main is still in progress at the edge of the gi¬ 
ant molecular shell of the complex. Using molecular emission maps, we showed that the 
two main gas clumps associated with these clusters (which together contain more than 
4.5x10^ Mq) show strong velocity gradients and gas outflow activity, spatially coincident 
with clusters and sub-structur es. This is in good agreement with the convergent flow scenario 
flRivera-Ingraham et al.ll2013l) that contributes to a large massive star formation efficiency 
in W3. Our analysis suggest that both the progression of formation following the primordial 
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structure of the original Giant Molecular Cloud, and the triggered formation caused by the 
compression and dispersal of gas (that is, within the cloud) are both important in W3. 

The W3 Complex, presents a rather complicated layout. The principal clusters overlap 
partially along the line of sight and they may still interact with each other. Moreover, the 
clusters present a signihcant sub-structure. If the sub-structures are bonahde sub-clusters, 
they pose an interesting problem for future higher quality data. Particularly, space based 
mid-infrared imaging with high resolution and deeper photometric limits (e.g. the James 
Webb Space Telescope) will provide more complete samples of the sub-cluster population. 
Also, spectroscopic data for individual young candidate members, may rehne age and age 
spread estimates. 
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A. Near-IR Observations 

In Table |3] we list information corresponding to our near-IR Calar Alto 3.5m OMEGA 
2000 observations. 


Table 3. Near-Infrared Observations of W3 Fields 


Field ID 

Date Obs. 

Center Coords. 

J2000 

Eilter 

Seeing 

[(")] 

BD Peak” 
[mag] 

OMEGA 2000-CAHA 3.5m OBSERVATIONS 

CAHA-W301 

2008-11-08 

36.580583 

62.135000 

J 

1.00 

21.00 

CAHA-W301 

2008-11-08 

36.580625 

62.135000 

H 

1.05 

20.00 

CAHA-W301 

2008-11-07 

36.580625 

62.135000 

K 

1.07 

19.50 

CAHA-W302 

2008-11-09 

37.100042 

62.108056 

J 

1.09 

20.50 

CAHA-W302 

2008-11-08 

37.100042 

62.108056 

H 

0.98 

20.25 

CAHA-W302 

2008-11-09 

37.100042 

62.108056 

K 

1.00 

19.25 

CAHA-W303 

2008-11-09 

36.659000 

61.985833 

J 

1.08 

20.75 

CAHA-W303 

2008-11-09 

36.659000 

61.985833 

H 

1.07 

19.75 

CAHA-W303 

2008-11-09 

36.658958 

61.985833 

K 

1.16 

19.25 

CAHA-W304 

2008-11-09 

37.056875 

61.904167 

J 

1.04 

20.75 

CAHA-W304 

2008-11-09 

36.580583 

62.135000 

H 

0.96 

20.50 

CAHA-W304 

2008-11-09 

37.056875 

61.904167 

K 

0.91 

19.50 

CAHA-W3CF’’ 

2008-11-09 

33.924208 

60.545000 

J 

1.07 

21.00 

CAHA-W3CF 

2008-11-09 

33.924208 

60.545000 

H 

1.14 

19.75 

CAHA-W3CF 

2008-11-09 

33.924250 

60.545000 

K 

1.06 

19.50 


’^Brightness Distribution Peak. Expresses turnover point of observed magnitude 
distribution 


’’Control Field 







